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Chiral d-wave superfluid is a preliminary example of topological matters that intrin¬ 
sically encodes interaction effects. It exhibits fascinating properties including a finite 
Chern number in the bulk and topologically protected edge states, which have been 
invoking physicists for decades. However, unlike s-wave superfluids prevalent in na¬ 
ture, its existence requires a strong interaction in the d-wave channel, a criterion that 
is difficult to access in ordinary systems. So far, such an unconventional superfluid 
has not been discovered in experiments. Here, we present a new principle for creating 
a two-dimensional chiral d-wave superfluid using periodically driven lattices. Due to 
an imprinted two-dimensional pseudospin-orbit coupling, where the sublattice index 
serves as the pseudospin, s-wave interaction between two hyperflne spin states natu¬ 
rally creates a chiral d-wave superfluid. This scheme also allows physicists to study 
the phase transition between the topologically distinct s- and d-wave superfluids by 
controlling the driving field or the particle density. 


Superfluidity takes center stage in modern condensed 
matter physics [TJ [2]. Whereas a large number of elements 
on the periodic table become superfluids or superconduc¬ 
tors at low temperatures, most of them are conventional 
s-wave ones. Though superfluidity could also occur in 
high-partial-wave channels, such as p-wave superfluids in 
Helium and d-wave superconductors in cuprates 13 SI, 
it is very rare that a high-partial-wave superfluidity is 
accompanied by the time-reversal symmetry breaking. 
Such chiral superfluids, for instance, the Px + ipy and 
dx 2 -y 2 + idxy superfluids [5HlQ]. give rise to exotic phe¬ 
nomena in an interacting many-body system. It has been 
realized that one need go beyond the conventional sym¬ 
metry breaking paradigm to describe these chiral super¬ 
fluids, since they are distinguished from their counter¬ 
parts that respect the time-reversal symmetry, such as 
Px and dx 2 -y 2 , by a finite Chern number, a topological 
property in the bulk, and chiral edge states protected 
by topology in a finite system. There have been experi¬ 
mental evidences for the Px + ipy superconductor, which 
exhibits a Chern number C = 1, to exist in Sr 2 Ru 04 
El- For the chiral dx2-y2 pidxy superfluid with a Chern 
number (7 = 2, there has been no experimental observa¬ 
tion yet. 

To realize a chiral dx‘ 2 -y 2 + idxy superfluid, theo¬ 
retical studies have proposed using doped graphene 
for enhancing the d-wave interaction na, or magnetic 
impurities for mixing dxy with dx‘ 2 -y 2 in cuprates [9}[T3]. 
These proposals requiring sophisticated doping of solid 
materials have not been realized in experiments so 
far. In ultracold atoms, though the interaction can 
be tuned by Feshbach resonance, the strong atom loss 
near high-partial-wave resonance makes it challenging to 
explore chiral superfluids with a large d-wave scattering 
length [m [T5] . Here, we propose a new scheme that 
simply requires 5-wave interaction. The idea is to 
engineer a special band structure in a lattice by breaking 


the time-reversal symmetry via a periodic driving 
field. Such a driving field imprints a two-dimensional 
pseudospin-orbit coupling in the single-particle band 
structure ^ kyUix + kxfUy, where k = (kx^ky) is the 
momentum, and the sublattice index plays the role of a 
pseudospin-1/2 m. Uploading two hyperflne spin states 
a =t, i that experience the same driving field to this lat¬ 
tice, the s-wave interaction between these hyperflne spin 
states inevitably induces a synthetic d-wave one, and the 
relative strengths of the interactions in the s- and d-wave 
channels can be well controlled by the driving field and 
other microscopic parameters. Since the time-reversal 
symmetry is readily broken in the single-particle level, 
when the synthetic d-wave interaction is dominant, a 
chiral d-wave superfiuid naturally emerges. 

Single-particle band structure We consider a two- 
dimensional lattice subject to a circularly polarized time- 
periodic gauge field A(r) = {Asinujr,—A cos ujt). The 
Hamiltonian in the Floquet framework can be written as 

Ho - ihdr = ^ ’’ + V{v) - ihdr, (1) 

where r = {x,y), and U(r) is the lattice potential. This 
Hamiltonian applies to both electrons in solids and ultra¬ 
cold atoms with synthetic gauge fields m- By a simple 
transformation x ^ x-\- f cos ujt, y ^ y-\- f sin ujt, t ^ r, 
one sees that equation 0 is equivalent to 

p2 

Hq — ihdr = + U(x -h / cos UJT, y f sin cor) — ihdr , 

( 2 ) 

where / = A/{Mu), except for a trivial position inde¬ 
pendent constant. Equation 0 describes a shaken op¬ 
tical lattice that has attracted a lot of interest recently 
in ultracold atom physics [17H2Q]. Our scheme there¬ 
fore works also for shaken optical lattices. Whereas our 
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FIG. 1: Periodically driven checkerboard lattice, (a) 

Applying a circularly polarized light (red helix) or a circular 
shaking (red circle) to the checkerboard lattice, the tunneling 
from a A sublattice site to its four nearest B sublattice sites 
becomes dependent on the direction, i.e., ~ diQ and ~ zbzQ 
along the x and the y directions respectively. Q is the strength 
of the photon-assisted inter-sublattice tunneling, (b) Red and 
blue clouds represent the Wannier wave functions in the A and 
B sublattice sites respectively. A is the energy offset of the 
two sublattices, and 5 — — hw is the one-photon detuning. 


scheme is general, to concrete the discussion, we focus on 
a checkerboard lattice composed of A and B sublattices, 

y(r) = -Vo (cos^(^) + cos^(^) + 2scos(^)cos(^)) , 

(3) 

where Vq is the lattice depth, d is the lattice spacing, and 
s is a tunable constant, as shown in figure 1(a). A strong 
lattice potential Vz{z) is applied along the ^ direction for 
creating a two-dimensional system. Whereas a checker¬ 
board lattice is common in solids, such an optical lattice 
can also be created for ultracold atoms [ 2 TJI 22 ]. The band 
structure of both the static and the dynamically driven 
lattices can be solved exactly using the plane-wave ex¬ 
pansion and the Bloch-Floquet theory respectively. For 
the static lattice, the lowest two bands are well separated 
from other ones for large Vb- The main contributions 
to these two bands come from A and B sublattices re¬ 
spectively, since a finite s produces an energy offset A 
between the nearest neighbor sites and suppresses the 
inter-sublattice tunneling, as shown in figure 1(b). How¬ 
ever, once such a lattice is periodically driven, the atomic 
orbital in A(B) sublattice sites could absorb (emit) a pho¬ 
ton and overcome the energy mismatch. Such a photon- 
assisted tunneling produced by the A(r) • P term hy¬ 
bridizes the lowest two bands if uj is close to their band 
gap, and the Hamiltonian in the rotating wave approxi¬ 
mation is written as 




(4) 


where e^k = 4dA^oskxdcoskyd + 2t'j^{cos2kxd + 
cos 2kyd) — S, €Bk = 4tB cos kxd cos kyd + 21'^ (cos 2kxd + 


cos2kyd), = 2fl{ismkxd — smkyd)e~^^^^ and 

(pk = diTg{{ismkxd — smkyd)e~'^^^^}. Ia (^b) and 

characterize the nearest- and next-nearest neigh¬ 
bors tunnelings in the A(B) sublattice, Q is the am¬ 
plitude of the photon-assisted intra-sublattice tunneling, 
and S = A — huj is the one-photon detuning. Though 
a finite s doubles the unit cell of the square lattice with 
5 = 0, we have defined kx and ky^ as shown in figure 2, to 
be consistent with the expression of the lattice potential 
shown in equation 

Equation 0 has a clear physical interpretation in the 
real space. As shown in figure 1(a), the inter-sublattice 
tunneling carries a nontrivial phase factor in the tight- 
binding model. For the tunneling from a A sublattice site 
to its four nearest neighboring B sublattice sites, the tun¬ 
neling is proportional to l,i,—1,—i respectively. These 
complex tunnelings come from the circularly driving field 
that break the time-reversal symmetry. The relevant 
basis wave functions are lTA(r)e*^^ and HB(r), in the 
language of Floquet theory, where ITA(r) and Wb{'^^) 
are the Wannier functions of the A and B sublattice re¬ 
spectively. Physically, the phase factor corresponds 
to absorbing one photon in the A sublattice. The tun¬ 
neling from A to B sublattice sites along the x and —x 
directions, which are referred to as are proportional 

to dr f drWB(r A dx) sm(u;T)Px(WA(r)e^^^), where 
X is the unit vector. Since both Wa{t^) and Wb{t^) 
have even parity, one naturally sees that t±x ^ As 
for the inter-sublattice tunneling along the y direction, 

t±y ~ - fo^ dr f drWB(r±dy) cos(ujT)Py(WA(r)e^^'^) ~ 

±if}. Fourier transform of the tight binding model to 
the momentum space leads to the off-diagonal terms 
in equation 0 with a momentum-dependent phase 
factor (pk- Near the F point, (fk ^ cirg{—ky + ikx}. 
This corresponds to a two-dimensional pseudospin-orbit 
coupling ^ kyirix + kxTriy^ where the sublattice index 
represents the pseudospin m. It is known that, if one 
uses the Raman scheme, phase-locked multiple lasers 
are required for realizing a multi-dimensional spin-orbit 
coupling. Here, circularly driving the optical lattice itself 
naturally provides such an important pseudospin-orbit 
coupling. 


Chiral d-wave superfluid induced by s-wave in¬ 
teraction Whereas the Hamiltonian in equation Q may 
lead to quantum anomalous Hall effect in the strong hy¬ 
bridization limit 4(^(4 — — ptA — tB\ < d < 4(f^ — 

t'^) 4\tA — we focus on the trivial single-particle 

band structure with Chern number C = 0 in this Article. 
Strikingly, interaction completely changes the topological 
properties of the system, and gives rise to a chiral d-wave 
superfiuid with a Chern number (7 = 2. To see this fact, 
we consider two species of particles interacting through 
5 -wave interaction, which are labeled as t and We first 
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FIG. 2: cZ-wave interaction induced by the phase wind¬ 
ing in the single particle band structure, (a) The purple 
arrows represent the pseudospin m formed by the A and B 
sublattices, and the blue arrows represent the projection of 
this pseudospin on the plane. At both the F and M points 
in the Brillouin zone, the pseudospin is perpendicular to the 
plane due to the absence of the band hybridization. The red 
and green spheres represent fermions with hyperhne spin up 
and spin down respectively. The dashed circle represents the 
Fermi surface and the shadowed region represents occupied 
states. Due to the phase winding in the single-particle spec¬ 
trum, the scattering acquires a d-wave component that de¬ 
pends on the angle between the momentum k and kb (b) 
Top view of (a). The projection of the pseudospin rotates 27r 
on the Fermi surface enclosing the F point, (c) 13 ^/in the 
Brillouin zone. /3k = 0 at the F and M points. 

diagonalize the single-particle Hamiltonian using, 

where a =t, is the eigenstate of the 

matrix in equation W. Both o^k and f3\^ are real, a^. 
and are the creation operators for the unhybridized 
bands. Due to the phase winding induced by the 
driving field, we refer the eigenstates as to helicity bands. 
Equation (§ avoids the ambiguity of the definition of the 
phase factor (pk at the F and M points, where the off- 
diagonal term in equation 0 vanishes and ((ak,/5k) = 
(1,0). In other words, there is no band hybridization at 
these two points in the Brillouin zone (BZ). If one treats 
the sublattice index as a pseudospin, such a pseudospin 
points along the same direction at both the F and M point 
in the same helicity band. Therefore, both bands are 
topologically trivial in the single-particle level with Chern 
number (7 = 0. However, in the presence of interaction, 
a topologically nontrivial superfluid rises as the ground 
state. 

Consider two hyperfine spin states a =t, i uploaded to 
the shaken lattice, the interaction Hamiltonian contain¬ 
ing only on-site attraction can be written as 

V = -Ua T] ’di-ffiii -Ub'^ ftitfiii, (6) 

ieA ieB 

where fiia- is the density operator for spin-cr particles 


at site i, and Ua,b > 0. Without loss of general¬ 
ity, we have chosen Ua Ub^ since the lattice po¬ 
tentials near the minima of the A and B sublattice 
sites are in general different. For ultracold atoms, 
the ratio UaIUb can be tuned in a broad range. 
First, for a fixed scattering length as{as < 0), Ur^ = 
J dz\ip{z)\^ f dr\4)ri{r)\^, where rj = A, B, 
and 07^ (r) are the Wannier wave functions of the lattice 
potential Vz{z) and E(r) respectively. As shown in figure 
1, A sublattice sites have lower energies than B sublattice 
ones for a finite 5 , and therefore tighter confinements. 
The Wannier wave function 0 a (i*) is thus more local¬ 
ized than 0 b(i*). By tuning s, the ratio Ua/Ub changes. 
Second, this ratio can be further tuned by introducing 
a position dependent as(r) using optical Feshbach reso¬ 
nance, since a standing wave naturally leads to a peri¬ 
odically modulating scattering length as(r) in the real 
space [23]. Alternatively, the current experimental ad¬ 
vancements on creating magnetic lattices allow one to 
produce a spatially periodic magnetic field with largely 
tunable lattice geometry, period and amplitude [ 2 ^ 1 ^. 
For a spatially dependent as(r), the on-site interaction 
then becomes Ur^ = —/ dz\(p{z)\^ f (iras(r)| 07 ^(r)|^. 
It is even possible to strongly suppress one of them if 
as(r) is tuned to zero at the lattice sites of one of the 
sublattices. Finally, there have been reports in a number 
of recent conferences that a well designed optical micro¬ 
scope generated by projecting a holographic mask [26| 
could create optical lattices with arbitrary lattice poten¬ 
tials, such as those going beyond traditional sinusoidal 
ones. Using this technique, 0 a (r) and 0 b (r) can be 
independent controlled and provide an extra degree of 
freedom to tune Ua/Ub- In this Article, we theoretically 
study how the ground state depends on Ua/Ub, and the 
driving field A(t). It is worth mentioning that all discus¬ 
sions here can be directly generalized to a related system, 
where the A and B sublattices are physically separated 
along the ^ direction so that ip{z) is different in A and B 
sublattices (Supplementary Note 1 ). 


We first consider a filling factor {fii^ + < 2 . 

In the weakly interacting regime, where the interaction 
strength is much smaller than the band gap between 
the two helicity bands, we project V to the lower he¬ 
licity band for constructing the low-energy effective the¬ 
ory. We concrete the discussions on positive detun¬ 
ing 6 > 4(t(^ — + 4|tA — ^bI, since all results can 

be straightforwardly generalized to negative detuning 
6 < — — 4|tA — ^bI- For a positive detuning, after 

absorbing a photon, A sublattice sites still have lower en¬ 
ergies than B sublattice ones. The main contribution to 
the lower helicity band comes from the A sublattice and 
the corresponding operator for the lower helicity band is 
d^.. The BCS type of projected interaction Hamiltonian 
is written as Ve/f = - 
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where 


T4k' = , 

(7) 

where N is the number of lattice sites. Equation 0 is a 
key result of our work. It reveals a significant effect of the 
band hybridization in the driven lattice. Though the bare 
interaction is purely a s-wave one, the phase winding en¬ 
coded in the two-dimensional pseudospin-orbit coupling 
^ kymx^kxfriy leads to an effective interaction that con¬ 
tains both the s- and d-wave components. This can be 
intuitively understood from the fact that both the a 
and a components see the same driving field. Since 
the single-particle spectrum of each hyperfine spin state 
exhibits a phase winding ^ the pairing scattering 

between a =t and a inevitably carries a d-wave com¬ 
ponent as shown in equation 0 . 

The relative strengths of interaction in s- and d-partial- 
wave channels can be well controlled by microscopic pa¬ 
rameters of the system, such as the ratio between the 
on-site interaction in A and B sublattices Ua/Ub- For 
the extreme cases, 1/^ = 0 or 1/^=0, the interaction 
is a purely s- or d-wave one respectively. It is expected 
that the chiral d-wave superfiuid shall emerge when Ub 
is dominant. Another important parameter is the ratio 
/^k/<^k, which is controlled by the single-particle spec¬ 
trum. At the T and M points where the band hybridiza¬ 
tion is absent, /3k = 0 and o^k = 1. The interaction is a 
purely s-wave one. Away from the T and M points, /3k 
becomes finite. The larger the driving amplitude A is, 
the larger the ratio /dk/o^k is. As shown later, the value 
of /dk/trk at the fermi energy is crucial to determine the 
topological properties of the superfiuid. 

Define the order parameters in the s- and d-partial- 
wave channels, 

As = —Ua <ak(<^-kt^kt)/A, 

(8) 

k 


the BCS mean field Hamiltonian can be written as 


Hm 



ek - 

A:k+A*i, 


Ask + A^k A 
/^ - ek y ’ 


(9) 


where ek and /i are the single-particle energy of the lower 
helicity band and the chemical potential respectively. 
Ask = AsO^k, and A^k = At small 

momenta, ^ iky)^. Thus A^ corre¬ 

sponds to the order parameter in the da, 2_^2 -\-idxy chan¬ 
nel. As and A^^ can be obtained self-consistently through 
the standard BCS approach (Supplementary Note 2), and 
the ground state wave function is. 


IG) = n ('“k + ^ke*®''4d-ki)|0)> (10) 

keBZ 



FIG. 3: Phase diagram, (a) The x and y axis represent 
the the ratio of the on-site interaction UaIUb and the fill¬ 
ing factor n respectively. Chiral d-wave, s-wave and the 
mixed partial-wave superfluids are represented by blue, red 
and green colors respectively. For a fixed n, the chiral d-wave 
superfiuid shows up at large enough Ub- The black curve 
represents the topological transition line characterizing the 
change of the Chern number from 2 to 0. For the numerical 
parameters we chose, it is very close to the boundary between 
the chiral d-wave and the s -h d superfiuid. (b) As and Ad in 
terms of the recoil energy Er = j{2MdP) as a function of 

UaIUb at half filling n = 1. Inset shows the region near the 
topological phase transition point highlighted by the dashed 
vertical line, (c) At a constant UaIUb-, different filling factors 
lead to different superfluids. On different fermi surfaces, the 
contribution from the A and B sublattice to the helicity band, 
which is characterized by ak/dk, is different, and different su¬ 
perfluids rise. (I), (II) and (HI) represent the fermi surfaces 
of the s-wave, the mixed, and the chiral d-wave superfluids 
respectively before turning on the interaction. 


where u\^ and v\^ are real, and dk = arglA^a^ -h 
Microscopically, A^ and A^ are con¬ 
trolled by UaIUb^ and the parameters of single-particle 
Hamiltonian, including both the detuning S and the 
strength of the driving field A. Two extreme cases are 
rather clear. When Ub = 0, A^ is always zero, and 
dk = 0. The ground state is a conventional s-wave super- 
fluid, |G)s = rikC'^k +'^k<3k>[^d^ki)l^)- contrast, when 
Ua = 0, As vanishes, and the superfiuid is dx‘ 2 -y 2 Eidxy 
superfiuid as expected, since the interaction in equation 
0 is a purely d-wave one. The ground state becomes 
\G)d = nk(^k + Using the 

pseudospin representation where and |0) are 

treated as the up and down component of a pseudospin 
Sz = ±1/2, Bloch sphere is covered twice by this pseu¬ 
dospin defined in the Brillouin zone, i.e., the Chern num¬ 
ber is 2, a characteristic feature of topological superfluids. 

For general values of Ua and Ub^ the phase diagram is 
shown in figure 3. For any given single-particle spec¬ 
trum with a fixed 6 and A, the chiral d-wave super- 
fluid is the ground state for large enough UbIUa- For a 
spatially dependent scattering length, we use the ansatz 
asir) = ao — ai cos(7rx/d) cos(7r^/d) to model it. For 
ao ^ 20nm, and ai ^ 16nm, which corresponds to 
UaIUb about 0.05 — 0.1, we find out that Ad can be 
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a few to ten nK. With further increasing Ub-, could 
be even larger. With the current advancements of cool¬ 
ing atoms [ 27 H 29 ] . it is promising to access such a tem¬ 
perature scale in experiments. For a uniform scattering 
length, one could also adjust Ub/Ua and access the chi¬ 
ral d-wave superfluid by adjusting the lattice potential 
itself via changing V and s (Supplementary Note 3). It 
is worth mentioning that, in an ordinary checkerboard or 
square lattice that respects the D/^ symmetry, it is diffi¬ 
cult to couple dx 2 _y 2 and dxy which belong to different 
irreducible representations of the D 4 group. In our sys¬ 
tem, the time-reversal symmetry is readily broken in the 
single-particle level so that such a constraint is absent. 
In the weakly interacting regime, particles prefer to stay 
in the lower helicity band, where the effective interaction 
inevitably includes a + idxy component, and the 

chiral d-wave superfluid naturally rises. 

With increasing UaIUb^ there is a first order transition 
from a pure chiral d-wave superfluid to a mixed s- and 
d-wave superfiuid in the mean field framework. Such a 
mixed superfiuid may be referred as to s-\-{dx2-y2 -\-idxy). 
Whereas both a pure s and a pure chiral d-wave super flu¬ 
ids are gapped, such a mixed state may be gapless by fine 
tuning of the parameters. When the gap closes, it cor¬ 
responds to a topological phase transition point defined 
by 

Asal,=Ad^l,, ek*=M (H) 


where k* is a momentum on the Fermi surface along the 
diagonal direction of BZ such that k* = 0. Along this 
direction, the off-diagonal term in the BCS Hamiltonian 
becomes real and could vanish, since ip = tt. 

When equation (11) is satisfied, both the diagonal and 
off-diagonal term of the BCS Hamiltonian in equation © 
vanish so that the excitation gap disappears. Across the 
transition point, and the Chern number of the ground 
state jumps from 2 to 0. On each side of the transi¬ 
tion point, the mixed state is topologically equivalent to 
either the chiral dx 2 -y 2 


idxy OT the s-wave superfiuid 


(Supplementary Note 4). It is worth mentioning that for 
the mixed s- and d-wave superfiuid, there is a two-fold 
degeneracy, since changing the sign of Ad in equation 
leads to the same ground state energy. For such a 
s — {dx‘ 2 -y 2 -h idxy) superfiuid, k* satisfies k* = 0 such 
that ip\^* -h ip-\^* = 0. 


Equation 0 tells one that, in addition to chang¬ 
ing Ua/Ub^ tuning the filling factor or the driving field 
also allows one to access the topological phase transition, 
since the amplitude of the driving field and the location 
of the Fermi surface determine the value of /dkp /<akF • 
shown in figure 3, for a Fermi surface very close to the F 
or the M point, /^kp/Q^kp is small, and the 5 - wave interac¬ 
tion is dominate so that the Chern number is 0. With a 
constant Ub/Ua^ enlarging the Fermi surface or increas¬ 
ing the driving amplitude A increases /3kp /<akp and leads 
to the transition to the ground state with (7 = 2. 


a b c 



fc^c?/(V27r) 


FIG. 4: Edge states, (a) When UaIUb is small, the topo¬ 
logically non-trivial ground state has two chiral edge states 
on each edge(red and green solid lines), (b) When UaIUb 
increases to a critical value, the gap in the excitation spec¬ 
trum vanishes. A topological phase transition occurs, and 
edge states disappear, (c) With further increasing UaIUb-, 
the gap reopens in the bulk and there is no edge state in the 
topologically trivial state. 


A characteristic feature of the chiral d-wave superfiuid 
is the presence of two chiral edge states, as shown in 
figure 4, which serve as unique signatures of the chiral 
d-wave superfluids in experiments. We have chosen the 
edge along the x' direction cutting A sublattices so that 
^'x ^ ^ ky is a good quantum number. Along the y' 

direction, we have applied the open boundary condition. 
In solid materials, such chiral edge states correspond to 
a quantized Hall conductivity (Jh = Ih. In ultra¬ 
cold atoms, such edge states could also be measured us¬ 
ing a variety of techniques [30], such as the expansion of 
the atomic cloud [31], and the generalization of the Bragg 
spectroscopy that is sensitive to angular momentum [32] . 
This allows one to directly visualize the edge states and 
identify the chiral d-wave superfiuid. Such edge states 
are topologically protected, since a small mixing with the 
5 -wave component in the s + (da, 2_^2 + idxy) superfiuid 
retains the Chern number (7 = 2. 

For a large filling factor {fii^ + > 2, the fermi 

surface shall form in the upper helicity band. One 
then needs to project the interaction to the upper he- 
licity band, and =-^ Ekk' '^kk'Ckfchv-k'^^k't, 
where 

= (UBalal + , 

( 12 ) 

similar to Veff for the lower helicity band, with inverted 
roles of (7 b and (7^, as well as those of a\^ and /3k• All 
previous discussions then apply. 

Discussions To reveal the underlying physics for cre¬ 
ating the chiral d-wave superfiuid, we have projected the 
interaction to one of the helicity band. Whereas this is 
well justified in the weak interacting limit, it is useful to 
investigate the interaction induced mixing between the 
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two helicity bands. We write the full Hamiltonian as 
-ff = Ek where a-^i, 6 _k 4 ,) 




CAk - Ai 




0 


— 1^ 

0 

A* 


Aa 

0 

/i - e^k 


0 

Ab 

/i - CBk 


/ 


(13) 

where e^k is the single-particle energy of = A, B 
sublattices, A^i = — X]k(^kt^-k^)/A and Ab = 

—[/b Z^k(^kt^-ki)/A characterize the on-site pairing 
amplitudes. For noninteracting systems, equation (13) 
describes two identical copies of equation (|^ for spin-up 
and spin-down particles. The presence of attraction leads 
to the particle-hole hybridization in both A and B sub¬ 
lattices. Using the basis of single-particle eigenstates as 
defined in equation (§ , the Hamiltonian can be rewritten 
as H = Ek^k-^/»^k, where d-ki, T’ 


and 


Mh = 


^ ^dk ^dk 0 ^k 

- e^k k 0 

0 k ^ck ^ck 

< 0 lljk M-eck/ 


o* 


(14) 


which includes both helicity bands where U^k = A^o^k + 
Uck = Asal + and 

z^k = <^k/^k(AAe“*^-*^ — Abc*^*^). For weak interactions 
that are much smaller than the gap between the two helic¬ 
ity bands, the block off-diagonal matrix can be ignored, 
and the upper and lower helicity bands can be investi¬ 
gated independently, as discussed before. It is easy to 
identify that = Ab, = Aa, Af = Aa, = Ab 
where the superscript c and d represent the results for 
two different helicity bands respectively. 

To include the effect of inter-helicity band coupling, we 
have numerically solved equation ( p!4| ), and verified that 
the mixing of the upper helicity band is exponentially 
small in the weak interacting limit. We have also verified 
that both the Chern number of ground state and the 
existence of the chiral edge states are not affected neither. 
This proves that the chiral d-wave superfluid discovered 
in this Article is stable against a weak inter-helicity band 
coupling. 

The study on topological matters is an important topic 
of current studies in both condensed matter and cold 
atom physics. Whereas a number of non-interacting 
topological matters, such as topological insulators, have 
been realized in experiments [331435] , the search for topo¬ 
logical superfluids, which naturally encode the interac¬ 
tion effects, has just started [36H39] . This Article has 
proposed a new principle for realizing a chiral d-wave 
superfluid without resorting to a strong d-wave interac¬ 
tion. With the induced d-wave interaction in the driven 
lattice, a bare 5-wave interaction is sufficient to create 
a chiral d-wave superfiuid. Since periodic driving fields 


allow one to engineer the Bloch-Floquet bands, it is ex¬ 
pected that applying such driving fields provides physi¬ 
cists unique opportunities to create a variety of synthetic 
high-partial-wave interactions. It is promising that even 
more exotic topological superfluids, as well as other in¬ 
teracting topological matters, shall be created. 

Our scheme also provides physicists an ideal platform 
for studying the interplay between the topology of single¬ 
particle band structures and interaction. Whereas signif¬ 
icant progresses have been made in understanding topo¬ 
logical matters in non-interacting systems, many ques¬ 
tions remain open when interaction is introduced. This 
article has focused on a case where the single-particle 
band structure is topologically trivial, with Chern num¬ 
ber (7 = 0. We have seen that the interplay between 
interaction and a special band structure that encodes a 
phase winding in the single-particle eigenstates gives rise 
to topologically non-trivial matters with Chern number 
(7 = 2 . Since the topology of the band structure in our 
system can be well controlled, we believe that this work 
will stimulate more studies on the interaction effect in a 
topologically nontrivial band structure of a periodically 
driven lattice. 
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Supplementary Figures 



FIG. 5: A body-centered cubic lattice with a tunable UaIUb • An additional lattice potential along the direction 
differentiates the lattices sites at the corner (red dots ) and the center (purple dot) of each cube, which form the A and B 
sublattices. In such a lattice, the lattice potential in different layers along the ^-direction is different, and the ratio UaIUb 
depends on the amplitude V^. For large enough such a three-dimensional lattice is decomposed to many independent 
two-dimensional systems, each of which is a checkerboard lattice with A and B sublattice shifted along the z direction. 



FIG. 6: Realize chiral d-wave superfluid with a uniform scattering length. Left, The lowest two bands (red and blue 
solid curve) of the static lattice with the parameter V = 6Er, s = 0.1. Red dashed line is the side band formed by the A band 
with one photon absorption. Middle, Floquet band structure after the lattice is periodically driven. Purple dashed line is the 
fermi energy. In this case, with a fixed scattering length as = 9.2nm, the ground state is a topologically non-trivial s + d wave 
superfluid. Right: the dependence of Ua/Ub on s with a fixed scattering length. 
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FIG. 7: Topological phase transition. (a,b,c) Small arrows characterize the strength and direction of the {Sx,Sy) 
component of the pseudospin, which represents the particle and hole states in the BCS wave function, in the s + d superfluid. 
Small filled and empty dot represent the location in the BZ where the pseudospin has only the Sz component with Sz = 1/2 
and Sz = “1/2 respectively. Along the black solid and purple dashed loop, the {Sx,Sy) component of the pseudospin has 
winding number 1 and —1. Blue solid line is the fermi surface and the shadowed region represents the occupied state. In 
(a). As is very small, two filled {Sz = 1/2) and two empty dots{Sz = “1/2) are present in BZ. The resultant superfluid is 
topological non-trivial(C = 2), since the Bloch sphere is covered twice by the pseudospin, (b) When As increases, the two 
filled dots move from F to M and touch the fermi surface, this is the critical point of topological phase transition, where the 
excitation gap vanishes at the positions of the red filled dots on the Fermi surface, (c) Further increasing As, no point in the 
BZ has a pseudospin state Sz = 1/2, i.e., the north pole of the Bloch sphere is not covered, and the superfluid is topological 
trivial(C = 0). (d,e,f) Corresponding excitation spectrum of the s + d wave superfluid, (e) The excitation becomes gapless on 
the fermi surface at the critical point of topological phase transition. 


Supplementary Note 1. An alternative scheme to realize a bipartite lattice with tunable Ub/Ua 
We consider a three-dimensional lattice, which may be created by four pairs of lasers, 

y(r) = -Vo {cos^{ — ) + cos^{-^)j -V{z)-sVocos{ — )cos{-^)-s \JVoV^ ^cos(-^) cos( —) + cos( —) cos( —) j , 

(15) 

where V{z) = Vz cos^(^) — cos(^) is a one-dimensional double well lattice along the 2 : direction. When V/ = 0, 
such a lattice is a body centered cubic lattice. A large enough separates the system into decoupled two-dimensional 
ones, each of which is composed of two square lattices separated along the 2 ; direction by a distance d, as shown in 
Fig.SI. View from the top, such a two-dimensional lattice is equivalent to a checkerboard lattice. Whereas all the 
results for both non-interacting and interacting systems are the same as the checkerboard lattice discussed in the 
main text, such a lattice has an advantage in practice. Since the A and B sublattices are readily separated along the 
2 ) direction, a scattering length as{z) that changes only along the 2 ; direction is sufficient, if one wants to use such a 
scheme to further change Ua/Ub^ 

Supplementary Note 2. Self-consistent solutions of the BCS Hamiltonian 

Define the order parameter: 


, = -UaJ 2 Ad = -Ub {d.^M/N 

k k 


A 


(16) 
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The mean-field interaction Hamiltonian is: 

kk' 

= ^ + (17) 

Define and ^k = ^k — where /i is the chemical potential. Using Bogoliubov 

transformation, one can get the self-consistent equation about the order parameter: 


A, 


Ad 



(18) 

(19) 


This set of equations is solved numerically for obtaining the results presented in the main text. 

Supplementary Note 3. Tuning UaIUb with a uniform scattering length 

In the main text, we have discussed how to tuning Ua/Ub using a spatially variant scattering length as(r) for 
accessing the chiral d-wave superfluid. Even with a uniform scattering length a^, the same goal can be achieved, since 
the phase transition from 5 to d-wave superconductor also depends on ydk/ok- By engineering the single-particle band 
structure so that ^dk is much larger than Ok at the Fermi surface, the d-wave interaction becomes dominant so that 
the chiral d-wave superconductor emerges. For the parameters shown in supplementary figure S2, we have found out 
that when the filling factor {fii^ -\- n^) = 1.6, A^ = O.OlFi^. 


Supplementary Note 4. Chern number of the (da^ 2_^2 idxy) superfluid 

For a pure d-wave superfiuid with A^ = 0, the mapping from the BZ to the Bloch sphere gives rise to a Chern 
number of 2. The T and M point correspond to the north and south pole, if one views the occupied and unoccupied 
states as spin-up and spin-down state of a pseudospin-1/2. This can be seen directly from equation (9) of the main 
text, since /3k = 0 at the T and M point. 

When As becomes finite, there will be four points in the BZ that lead to vanishing off-diagonal terms of the matrix 
in equation (9) of the main text for a given ratio Ag/Ad. For 5 ± (da, 2_^2 -h idxy) superfiuid, this corresponds to 
AgOf^ ~ AAc^/3^. As shown in supplementary figure 3, if none of these four points touches the Fermi surface, the 
spectrum remains fully gapped, since the matrix is always finite. Moreover, it is easy to see that the Bloch sphere is 
still covered twice, and the Chern number remains to be 2 . Under this situation, the mixed superfiuid is topologically 
equivalent to a chiral d-wave superfiuid. With increasing A^, two of these four points will touch the Fermi surface, 
and at these two points, the excitation spectrum becomes gapless. If A^ further increases, all these four points will be 
located either within or outside the Fermi surface, so that either the south or the north pole cannot be covered. This 
means that the Chern number becomes zero, and the mixed superfiuid is equivalent to a trivial 5 -wave superfiuid. 



